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ABSTRACT 

We report on a recent spectral line survey of the planetary nebula (PN) 
NGC 7027 using the Arizona Radio Observatory (ARO) 12 m telescope and 
the Heinrich Hertz Submillimeter Telescope (SMT) at millimeter wavelengths. 
The spectra covering the frequency ranges 71-111 GHz, 157-161 GHz, and 218- 
267 GHz were obtained with a typical sensitivity of rms < 8mK. A total of 67 
spectral lines are detected, 21 of which are identified with 8 molecular species, 32 
with recombination lines from hydrogen and helium, and 14 remains unidentified. 
As the widths of emission lines from CO, other neutral molecules, molecular ions, 
as well as recombination of H + and He + are found to be different with each other, 
the line strengths and profiles are used to investigate the physical conditions and 
chemical processes of the neutral envelope of NGC 7027. The column densities 
and fractional abundances relative to H2 of the observed molecular species are cal- 
culated and compared with predictions from chemical models. We found evidence 
for overabundance of N 2 H + and underabundance of CS and HNC in NGC 7027, 
suggesting that X-ray emission and shock wave may play an important role in 
the chemistry of the hot molecular envelope of the young PN. 

Subject headings: ISM: molecules — radio lines: ISM — ISM: abundance - 
planetary nebulae: individual (NGC 7027) 
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Introduction 



We now realize that the circumstellar envelopes of evolved stars are active sites of molec- 
ular synthesis. Millimeter-wave spectroscopy has identified over 60 molecular species in the 
envelopes of asymptotic giant branch (AGB) stars, their descendants planetary nebulae 
(PNs), and the transition objects betw e en the two phases , proto-planetary nebulae (PPNs) 
( jOlofssonl 119971 ; ICernicharo et al.l |2000| ; IZiurys et al.l 120071 ) . These gas-phase molecules are 
believ ed to be the precursors of complex org anic compounds that are found in PPNs and 
PNs flCernicharo et al.l boOlal lbl: 1 Kwokl 12004 ). Since the molecular synthesis takes place in 
the circumstellar envelope through either gas-phase or grain-surface reactions, the study of 
the changing chemical composition and molecular abundance between objects in consecutive 
phases of evolution will provide useful information on the chemical pathways of their forma- 
tion. Furthermore, since the circumstellar envelopes are expanding, the dynamical timescales 
impose a time limit on the reaction time scales. The typical dynamical timescales for AGB 
stars, PPNs, and PNs are 10 4 -10 5 yr, 10 3 yr, and 10 3 -10 4 yr, respectively. The observation 
of a specific molecular species in the PPN phase which is not seen in AGB stars implies 
that this molecule must be formed on timescales shorter than a few thousand years. The 
changing relative abundance of molecular species between different phases can also point to 
specific chemical pathways. Motivated by these considerations, we have performed spectral 
line surveys of several objects in the late stages of stellar evolution and try to provide a 
systematic and uniform study of circumstellar chemistry. 

The study of the circumstellar chemistry of PNs is particularly interesting and challeng- 
ing because of the additional dimension of photochemistry. Since part of the circumstellar 
envelope is ionized by ultraviolet radiation from the hot central star, molecules can be 
photoionized or photodissociated. The molecular composition of the circumstellar envelopes 
is therefore expected to be qualitatively different from AGB stars and PPNs. The young, 
carbon-rich PN NGC7027 is known to be the PN with the richest molecular content. Be- 
cause of its high-mass (0.7 M ) central star, the central star of NGC7027 is evolving very 
quickly, reaching an effective temperat ure of 219 000 K in a few hundred years since the onset 
of photoionization fjZhang et al.ll2005l ). The remnant of the AGB envelope is eroded by a 
nascent ionized region of 10". Imaging observations have shown that a thin photodissocia- 
tion region (PD R) separates the ionized r egion and the extended molecular envelope with 
a size > 70" (see iHasegawa fc Kwokll200ll . and the references therein). A large number of 
molecular emis s ion lines have b e en detected by infrared ( I R) and radio spect roscopy (e.g. 
Liu et all Il996l : ICox et all 1 19931 : Irlaseeawa fc Kwokl l200ll ; Irlerpin et all l2002h . Unambigu- 
ously detected molecular species in NGC7027 include H 2 , OH, CO, CO+, CH, CH+,HCN, 
HNC, HCO+, H 2 0, N 2 H+, CN, CS, C?H C3H,, SiS, and their isoto pomers. A system- 
atic molecular-line observations of PNs by iBachiller et al.l (119971 ) and iJosselin fc Bachiller 
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( 120031 ) has found that the column densities of molecules observed in NGC 7027 are generally 
larger than those in other PNs. Compar ed with other PNs, NGC 7027 shows a low HNC 
abundance and an abnormal high N 2 H + f lBachiller et al.lll997l ). lHasegawa fc Kwokl (120011 ) 
suggested that most of the observed molecular lines are emitted from the PDR of the ob- 
ject excep t for C O emission which arises from a much more extended region. According to 
Cox et al.l (119871). lines from the ri n g mol ecule C3H2 are emit t ed fro m dense regions. The 
observations by iThronson fc Ballyl (119861 ) and iDeguchi et al.l (119901 ) suggested that grain 
chemistry probably plays an important role for the formation of molecules in NGC 7027. 



The chemistry in the PDR of NGC 7027 has been modeled by lYan et al.l (119991 ). who 
suggested that the molecular abundance is significantly influenced by both photodissocia- 
tio n and shock chemistry. A detailed chemical model of the neutral envelope of NGC 7027 
by lHasegawa et al.l (120001 ) suggested that most of the detected molecules in NGC 7027 are 
formed over a short ti mescale after the obj ect evolved into the PN stage. This was supported 
by the observations of iHerpin et al.l (120021 ) with the Infrared Space Obse rvatory (ISO). How- 
ever, different conclusion was obtained by IJosselin fc Bachillerl (120031) who suggested that 



most of molecules can survive during the transition from AGB to PN. iRedman et al. I (120031 ) 
modeled the chemical evolution from AGB to PPN and PN considering the affect of dumpi- 
ness. By comparing model results and observations, they proposed that the clumps form at 
an early stage of PN and may be plentifully present within NGC 7027. 

Hydrogen and helium recombination lines from NGC 7027 are detectable at millime- 
ter wavelengths, allowing us to investigate the physical conditions and dynamics of the 
ionized regions. An advantage of the radio data is the reduced effect of dust extinc- 
tion. A few radio recombination lines from H56a-H110a have been detected in NGC 7027 



( IRoelfsema et al 



1991 



Chaisson fc Malkan 1976; Ershov fc Berulis 1989. and the references 



therein). lErshov fc Berulis I (119891 ) constructed a dynamical model to explain the observed 
radio recombination lines and found that the ex pansion velocity o f the i onized gas is in good 
agreement with that of the molecular envelope. IRoelfsema et al.l (Il99ll ) observed the H76a 
line observations and argued that the ionized gas is decelerated by the surrounding molecular 
material. 

In this paper, we report a systematic spectral survey of NGC 7027 using the ARO 12 m 
and the SMT 10 m telescopes □, covering the frequency ranges 71-111 GHz, 157-161 GHz, 
and 218-267 GHz. These data allow us to investigate the physical conditions and chemical 
compositions of the molecular envelope of this PN. This paper is organized as follows: Sect. 2 



1 The 12 m telescope and the Heinrich Hertz Submillimeter Telescope (SMT) is operated by the Arizona 
Radio Observatory (ARO), Steward Observatory, University of Arizona. 
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describes the observations and the procedures of data reduction; Sect. 3 presents the identi- 
fication and measurements of detected emission lines; in Sect. 4 the method of data analysis 
is described and the resulting excitation temperatures, column densities and abundances 
relative to H 2 are presented; in Sect. 5, we discuss the physical and chemical characteristics 
in the molecular envelope and ionized regions; and the main conclusions are summarized in 
Sect. 6. 



2. Observations and data reduction 

The spectral survey was carried out between 2005 April and 2006 September. The 
observations were made in beam switching mode with an azimuth beam throw of 2'. Pointing 
accuracy was checked every two hours. The on-source integration time was more than one 
hour for each frequency setting. 

The 71-111 GHz and 157-161 GHz spectra were obtained using the ARO 12 m tele- 
scope at Kitt Peak. The dual-channel SIS receivers were employed in single sideband mode, 
yielding a system temperature of 150-400 K. The image rejection ratio was 16-20 dB. The 
receiver back-ends were equipped with two 256-channel filter banks (FBs) with spectral res- 
olutions of 500 kHz and 1 MHz and a millimeter autocorrelator (MAC) with 3072 channels 
and 195 kHz resolution. All spectrometers were operated in series mode. Two orthogonal 
linear polarization modes were measured. The temperature scale, was calibrated and 
corrected for atmospheric attenuation, radiative loss, and rearward and forward scattering 
and spillover by the chopper-wheel method. The main beam brightness temperature was 
derived using Tr = where nf m is the corrected beam efficiency Over the frequency 

range 71-161 GHz, rf m is from 0.94 to 0.71, and the conversion factor is 32.5-29.1 Jy/K. 

The 218-267 GHz observations were carried out with the SMT 10 m telescopes at Mt. 
Graham, Arizona, using the dual-channel SIS receivers operated in single sideband dual 
polarization mode. The system temperature is typically 400-700 K. The spectrometers used 
were a 2048-channel acousto-optical spectrometer (AOS) with a spectral resolution of 500 kHz 
per channel and 1024-channel Forbes Filterbanks (FFBs) with a spectral resolution of 1 MHz 
per channel. All spectrometers were used simultaneously. The data were calibrated to the 
antenna temperature scale, T|, which was corrected for atmospheric attenuation. T\ is 
converted to main beam temperature by T R = T^/rj m b, where the main beam efficiency, 
r] m b, is ~ 0.7. The conversion factor is 35 Jy/K. 
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All spectra were reduced using the CLASS software package in GILDAS q After dis- 
carding the bad scans which are seriously affected by bandpass irregularities, we co-added 
the calibrated spectral data using the rms noise of each spectrum as weights. A low-order 
polynomial baseline defined by the line-free spectral regions was subtracted from each spec- 
trum. In order to improve the signal-to- noise, the ARO 12 m and SMT 10 m spectra were 
smoothed and rebinned to a spectral resolution of 1 MHz and 3 MHz, respectively, yielding 
a typical rms noise temperature of < 8mK in main beam temperature unit. Based on a 
few strong lines detected in different spectrometers and different epochs, we estimate that 
the calibration uncertainties introduced by instruments and atmospheric conditions amount 
to about 15%. 



3. The spectra 

The spectra of NGC7027 obtained with the ARO 12 m telescope and the SMT 10 m 
telescope are presented in Fig. [TJ and Fig. [2J respectively. The spectra at a full resolution 
are shown in Fig. [3] and Fig. HI which is only available in the electronic version of this 
manuscript. Note that some features in these figures are caused by bandpass irregularities. 
More molecular lines were detected in the 10 m spectra compared to the 12 m spectra as 
the ARO 12 m telescope has a larger beam size and thus bears larger beam dilution effect. 
Sections of the spectra containing spectral lines are expanded in Figs. (12 m) and M (10 m) 
in order to show the line profiles. Since each spectrum was simultaneously observed with 
two receiver back-ends, an emission line should be recorded by both spectrometers. This 
allows us to distinguish between real emission lines and artificial features. Contamination 
from image band was carefully checked. The feature at 221.463 GHz was found to be the 
CO (2-1) line from the image band. No other lines from the image band were apparent. A 
total of 71 distinct emission features were measured. 



The molecular lines are identified on the basis of the JPL c atalog (IPickett et al. 



1998] 



and the Cologne database for molecular spectroscopy (CDMS, iMuller et al.l 120011 . l20051 n. 
The suggested identifications are labeled in Figs. [Hand [2j Combining the single and blended 
features, we identified 21 molecular lines, belonging to 8 molecular species (CO, CN, C2H, 
HCN, HCO+, HCS+, HC 3 N, N 2 H+) with a total of 13 isotopomers. We also present 
'possible' detection of C3H2, HC 17 + , and HC 18 + . 14 lines remain unidentified. Several 



2 GILDAS is developed and distributed by the Observatorie de Grenoble and IRAM. 

3 http://spec.jpl. nasa.gov. 

4 littp: / /www. phi .uni- koeln.de /vorhersagen/ . 
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of these transitions are discovered for the first time in NGC 7027. A list of molecular lines 
detected in our spectra is presented in Table HJ The table also gives the rms noise, the main 
beam temperatures, and the integrated intensities (in units of K km/s) of the lines. 

The CO (2-1) and 13 CO (1-0, 2-1) lines have been ext ensively observed by vari- 
ous telescopes in the past (see, e.g. lHasegawa fc Kwok! l200ll ). The integrated intensi- 
ties of th e CO (2-1) and 13 CO (2- 1) lines presented in Table [1] are lower than those re- 



ported by lHasegawa fc Kwok! (120011 ) (see their Table 1). However, we note that the CO(2- 
l) / 13 CO(2-l) integrated int ensity ratio of 24.8 is in excellent agreement the value of 21.8 



in 



Hasegawa &: Kwok! (120011 ) . The lower antenna temperatures detected here are due to the 
larger beam dilution effect since the beam size of the SMT 10 m telescope is a facto r of_~ 2 



larger than the 15 m James Clarke Maxwell Telescope used by lHasegawa fe Kwok! (120011 ). 
The C 18 (2-1) and C 17 (2-1) lines were clearly detected in our spectra. To our knowledge, 
it is the first time that these two transitions are detected in NGC 7027. 

Three strong CN (2-1) hyperfine group s are detected in our o bserv ation s. These lines 
have p reviously been seen by lCox et al.l (119931 ). I Josselin fc Bachillerl (120031 ) . and lBachiller et al 
(119971 ). The emission feature at 227.2 GHz is only marginally above our limit of detection 



and may be a combination of the CN (N = 2 — 1, J = 5/2 — 1/2) and the C3H2 (432-821) 
transitions. Other previous detections of CN i n NGC 7027 include the J = 3 — 2 and J = 1 — 



transitions at 340 and 114 GHz, respectively (lHasegawa fc Kw ok 200ll; iBachiller et al.lll997 
Thronson fc Ballv!ll986l : ICox et al.lfl993l : I Josselin fc Bachillerl booj l. 



The C 2 H (3-2) transitions in NGC 7027, first discovered bv lHasegawa & Kwokl (l200ll ). 
are confirmed in our observations. All the C2H (3-2) lines show a double-peaked profile. 
The C3H2 (22o - 2n) at 18.3 GHz a n d the C3H2 (lio - loi) a t 21.6 GHz transitions were first 
detected in the PN by ICox et al.l (119871 ). We attempted to search for two new C 3 H 2 
transitions (6 2 5-5i 4 and 5 33 -4 2 2) at higher frequencies. The two lines are extremely faint, 
with intensities between 2-3<r noise level. Hence, our detection of C3H2 should be treated 
with some caution. 

The well studied HCN (1-0, 3-2) and HCO+ (3-2) transitions are also detected. We 
also detected their isotopic transitions H 13 CN (3-2) and H 13 CO + (3-2). However, the 
transitions H 13 CN (1-0) at 86.340 GHz and H 13 CO+ (1-0) at 86.754 GHz are overwhelmed 
by noise. The transitions HC 17 + (1-0) and HC 18 + (1-0) are only marginally above 2a 
noise level, and thus have an ambiguous detection. A narrow feature at 256.028 GHz has 
be en discovered and identifie d with the HCS + (6-5) transition, confirming the observations 
by 
by 



Cox et al 



Jl993h . 



Hasegawa fc Kwokl (1200 ll ). We also detect the N 2 H + (1-0) line which was first reported 
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An unsuccessful sear ch for the HC3N (11-10) transition at 100.076 GHz was made by 
Thronson &: Bally ( 19861 ) Fl In our observations, several HC3N lines are identified in the 
spectra from both the ARO 12 m telescope and SMT 10 m telescope. However, we do 
not detect the transitions HC 3 N (10-9) at 90.979 GHz and HC 3 N (12-11) at 109.174 GHz, 
making our identification arguable. If confirmed, this is the first detection of this species in 
a PN. 

The unidentified lines are listed in Table [2J Most of them are very faint, and are only 
marginally detected. S ome of these lines have been detected in the star formation regions Sgr 
B2(N) and Sgr B2(M) (jTurnerlll989l ; iNummelin et al.lll998l ). We also compared our spectra 
with a recent molecular line survey of the PPN CRL618 by IPardo et al.l (120071 ). These U- 
lines are not detected in their spectra. Thus we infer that some of these lines probably arise 
from species which are formed in very strong UV radiation. 

Table [3] lists a number of molecular species which are not detected by our observations 
of NGC 7027. For a given species, the strongest transition in our frequency ranges is given 
in the table aiming to reach the tightest constrain on its column density. The upper limits 
of the integrated intensities can be obtained by 3aAv/y/n, where a is the rms noise, and n 
is the number of channels within the line width Av. 

A number of recombination lines, including 24 hydrogen lines and 8 helium lines, are 
detected in our observations. Previously, most of the recombination lines were observed in 
the low- frequency bands. The id entifi cations of recombin ation lines are based on the calcula- 
tions by iLilley &: Palmer I (119681 ) and lTowle et al.l (119961 ) . No heavy-element recombination 
line is visible in our spectra as they are under our detection limit. The identifications and 
measurements of the recombination lines are presented in Table HI All recombination lines 
were measured using Gaussian line profile fitting. Some of the lines, such as H36/3, are prob- 
ably blended with unknown features and show a relatively broad profile. The H6O7 line at 
84.914 GHz The H6O7 line and The He42a line at 85.723 GHz fall within the spectral regions 
which have a low signal-to- noise (see Fig. [1]), and thus their measurements are uncertain. 



5 Note that in Table 1 of 
as "110.0 GHz" by mistake. 



Thronson &; Ballv 



(1986), the frequency of the HC3N (11-10) transition is given 
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4. Analysis 

4.1. Excitation temperatures and column densities 

In order to obtain molecular excitation temperatures and column densities, the effect of 
beam dilution is corrected under the assumption that both the source brightness distribution 
and the antenna beam have a Gaussian profile. The source brightness temperature (which is 
equivalent to source intensity) was obtained by T s = Tr(91 + 6 2 s )/0 2 s , where the main-beam 
brightness temperature Tr has the physical meaning of average intensity, 9^ is the antenna 
full beam at half-power and 9 S is the source diameter. The beam sizes of the ARO 12 m 
telescope and the SMT 10 m telescope are 86"-38" and 32"-28" over the frequency ranges 
71-161 GHz and 218-268 GHz, respectively. Different molecular species in NGC7027 might 
have different distribution and 9 S thus should be separately estim ated in terms of map of each 



species. Based on the observations of lHasegawa fc Kwokl (]200ll ). we have assumed 6 S = 13'! 2 



for all the species except CO. High-resolution mapping observatio ns of CO in NGC 7 027 



suggest that the CO emission originates in a more extended region (IMasson et ajj [19851 ), so 
we adopted 9 S = 40" for CO and its isotopomers. 

The standard rot at ion- diagram method was applied to determine molecular excitation 
temperatures and column densities. Under the assumption that i) the lines are optically 
thin, ii) the level populations are in LTE and can be described by a Boltzmann tempera- 
ture, referred to as the "excitation temperature" (T ex ), and iii) T ex is much larger than the 
background temperature (T bg ), we have the relation 

, N u , 3k f T s dv N E u 

In — = In = In - — - , 1 

g u 8tHz/6/^ Q{T ex ) kT ex 

where N u , g u , and E u is the population, degeneracy, and excitation energy of the upper level, 
J T s dv is the integrated intensity of the source, S the line strength, /i the dipole moment, 
v the line frequency, Q the rotation partition function, and N the column density. If many 
transitions for an individual molecule are observed, a straight line can be fitted by plotting 
the integrated intensity versus upper level energy. Then T ex and the total column density 
can be deduced from the slope (— 1/T ex ) and intercept [ln(A/Q)] of this line. 

Albeit with some problems, we probably detect five HC3N transitions. The excitation 
energies of the upper levels cover a large range with 10 K< E/k < 180 K, allowing us to 
determine excitation temperature and column density. Fig. [7] shows the rotation diagram 
for HC3N. A linear least-squares fit to all the transitions produces N = 4.39 x 10 13 cm~ 2 . 
However, the data in the rotation diagram s uggest that the low E„ transitions seem to arise 



from a colder region. The results presented in lHasegawa et al.l (120001 ) suggest that the neutral 
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envelope of NGC 7027 consists of a geometrically thin and dense shell and an extended outer 
stellar wind region, with the temperature in the dense shell being more than a factor of 10 
higher than that in the wind region. In this picture, our HC3N high E u transitions would 
originate dominantly in the hot, dense shell, whereas the cold stellar wind region contributes 
significantly to the excitation of low E u lines. The stratified temperature structure may cause 
some problems in the calculations of column d e nsity. To get more precise result s , deta iled 
modelling is required, as done by iPardo et al.l (12004 ) and iPardo fc Cernicharo I (120071 ) for 
CRL 618. However, this is hampered by the weakness of most molecular lines in NGC 7027. 

No rotation diagram is presented for other molecular species because either only one 
line for the species was detected, or the upper levels of the observed transitions have similar 
excitation energies. We thus adopted a constant excitation temperature of T ex = 34.6 K, 
as deduced from the HC3N rotation diagram, and applied Equ. ([T]) to estimate the column 
densities of CN, C 2 H, C 3 H 2 , HCS+, N 2 H+, and HCO+. We derived a column density of 
iV = 1.02 x 10 14 cm -2 for CN, which is consist e nt wi th th e values of 3.5 x 10 1 4 cm" 2 and 



6.1 x 10 13 cm~ 2 deduced by lHasegawa fc Kwokl (1200 ll ) and IJosselin fc Bachillerl (120031 ) . re- 
spectively. Our analysis yields N = 3.49 x 10 14 cm ~ 2 for C 2 H, in good agreement with 
the value of 2.4 x 10 14 cm -2 given by lHasegawa et al.l (J2000J). Fro m the obser v ations of the 
C 3 H 2 (lio-loi) at 18.3 GHz C 3 H 2 (2 2 o-2 1:1 ) at 21.6 GHz transitions, ICox et all Jl987h found 



iV(C3H 2 )~ 3 x 10 cm" , which is lower than our value (1.11 x 10 cm" ) by a factor of 
about 36. Given large measurement errors, our value may be unreliable. For HCS + , we 
obtain a value of N = 4.96 x 10 12 cm -2 , which should be treated with some caution as it 
was derived from an extremely faint emission line. From the N 2 H + (l-0) line at 93.2 GHz, 
we derive iV(N 2 H + )= 2.43 x 10 13 cur 2 . The column densities of HCS + and N 2 H + have 
not been reported previo usly. Our determination of iV(HCO + ) is higher than that of 
Hasegawa fc Kwokl (120011 ) by a factor of 5. Considering the uncertainties in calibrations, 
such difference is not unreasonable. 

The strong HCN and CO transitions are likely to be optically thick and therefore the 
rotation diagram analysis is not applicable. For these lines we used the following equations 
to deduce their excitation temperatures and column densities. 

The equation of radiative transfer is 



T s =[I(T ex )-I{T bg )](l-e 



(2) 



where I(T) = hu[k(e hu ^ kT — 1)] _1 with h the Planck constant and k the Boltzmann constant, 
t is the optical depth at the center of a line, and Tb g is the cosmic background temperature 
(2.7K). For the optically thick lines (r ^> 1), T ex can be obtained from Equ. ([2]) and is given 

by 
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T 

J. p, 



fcln 



hv 



hv 



k[T s +I(T bg )} 



1 



(3) 



Assuming that the isotopomer has the same excitation temperature as the main line, 
we can obtain its optical depth at the line center from Equ. (0) 



n = In 



I(T ex )-I(T b9 )\ ■ 



(4) 



The column density of the isotopomer can be derived from 



Ni 



3hAvnQ(T ex )eMw". 



87r 3 / u 2 S'[exp( 



hv 

^ L ex 



1] 



(5) 



where Av is the line width. 



Assuming that lines of the isotopomer is optically thin, the T s ratio of the main line 
and the isotopic line can be derived using Equ. (E]): 



1 — exp(— t„ 



(6) 



Once Tj is obtained from Equ. (jl]), the optical depth of the main line, r m , can be obtained 
from the above equation. 



Finally, we can determine the column density of the main species using 

N r 



(7) 



Applying the above equations to HCN and CO, we obtained their excitation tempera- 
tures and column densities. The results are summarized in Table El O ur determination of 
the co lumn density of HCN is a factor of 16 higher than that derived by lHasegawa fc Kwok 
(120011 ) . In our calculations, a large error may be caused by the very faint H 13 CN line. If 
applying Eq u. (HI), we derive a l o wer li mit of iV(HCN)= 4.2 x 10 13 cm~ 2 , in agreement with 
the result of lHasegawa fc Kwotd (]200ll ). 



For compar i son, Table [5] also lists the column densities predicted by the model of 
Hasegawa et all (120001 ) . The column densities deduced by current work are in of the same 
approximate orders of magnitude with the model predictions except for CO and N^.H " |". Th e 
model overestimation of the column density of CO is addressed by lHasegawa et al.l (120001 ) . 
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and an improved model (IHasegawa fc Kwokll200ll ) gives a better fit to the observed column 
density of CO. The model, however, fails to explain the extremely high column density of 
N 2 H + . A detailed discussion will be presented in the next section. 

Assuming a line width of 30 km s -1 , we estimate the upper limits of the column den- 
sities for several non-detected molecular species, and these values are presented in Table [3j 
Since these upper limits may still serve as useful constraints of the chemi cal processes in 



NGC7027, these upper limits are compared with the model predictions of IHasegawa et al. 



( 120001 ) . Table [3] shows that the observed upper limits are in reasonable agreement with the 
model predictions except for CS, which obs erved column densi t y is l ower than the model 
prediction. Even in the improved model of IHasegawa fc Kwokl (120011 ). the theoretical col- 
umn density of iV(CS) = 8.4 x 12cm~ 2 is still higher than the upper limit derived in the 
present study. 



4.2. Abundances calculations 

In order to calculate the fractional abundance of a molecular species relative to H2, we 
need to know the number density of the species and that of H2, which can be derived by 
solving the coupled radiative transfer and statistical equilibrium equations. A full treatment 
of the line excitation requires knowledge of the geometrical structure of the envelope, the rates 
for collisional excitation of the molecules, and the effects of dust emission and absorption, 
which however is poorly known, and is beyond the scope of this paper. Instead, assuming that 
all molecular emission originate in a spherical sh ell, we ca l culate d the molecular abundances 
respect to H 2 using the expression suggested by Olofsson ( 19961 ). 



/ x = 1.7x10" 



_ 28 v e 9 b D Q(T ex )u 2 ul e E '/ kT - J T R dv 



M H2 9uA ul J^e-^ 2 dx' 

where the integrated intensity J T R dv is given in Kkms -1 , the full half power beam width 
6b is in arc sec, v e is the expansion velocity given in kms -1 , D is the distance in pc, Mh 2 is 
the mass loss rate in M yr _1 , v u \ the line frequency in GHz, g u the statistical weight of the 
upper level, A u i the Einstein coefficient for the transition, E\ the energy of the lower level, 
and Xi >e = Ri >e / (9f,D) with Ri and R e the inner radius and outer radius of the shell. The 
equation is deduced under the assumption that the lines are optically thin, T ex is uniform 
throughout the shell, the molecular density follows an r~ 2 law, and the shell is formed by a 
constant mass loss rate and has a constant expansion velocity. One should bear in mind that 
the spatial variations in the fractional abundances and excitations may be a large source of 
error. For the calculations, we adopted the excitati on temperatu res deduced in section 4.1 



(see Table E) and a distance 880 pc to NGC7027 (iMassonl Il989h . Following the model of 
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Hasegawa et al.l (120001 ). we assumed Ri = 0.017 pc and R e = 0.14pc. The expansion velocity 



was determined from the widths of a few strong lines, which suggest v e ~ 30kms 1 . 

Assuming that the CO (2-1) emission arises from the stellar-wind regi on, we can es- 



timate the mass loss rate from the line strength using the formula given by IWinters et al. 



d2002|) 

M — K 7 v 1 n~ 20 . 

8{J)fco 



M = 5.7xlO- 2O ^ o 2 f M yr- 1 , (9) 



where the correction factor S(J) = 0.5 for the CO (2-1) line, the distance D = 880 pc, the 
abundance of CO relative to H 2 , /co> is assumed to be 1 x 10~ 3 , and the expansion velocity 
v e = 30kms _1 . The derived mass loss rate of 1.1 x lO _4 M yr _1 is in good acc o rd wi th the 



value of 1.4 x 10 4 M Q yr 1 derived from the model calculations by ISopka et al.l (119891 ). The 



error in the determination of M due to the distance uncertainty is 34%. 



Using Equ. (jSJ), we derived fx for all the observed molecular species and the results are 
given in Table The calculated abundance ratios between different species are not rigorously 
equal to the ratios of column densities. This is due to two reasons. One is the different 
assumptions made to deconvolve the brightness temperature distribution of the source from 
the antenna beam. For the calculations of column densities, we have simply assumed that 
the source brightness distribution is Gaussian. Second, for the abundance calculations, we 
have assumed that the lines are optically thin. Consequently, when the emission is optically 
thick, fx given in the table should be treated as a lower limit. We can use the N(X)/fx 
ratio, which approximately gives the column density of H2, to estimate the error due to 
these assumptions. The last column of Table [5] gives the N(X)/ fx ratio deduced from each 
molecular species. These values are generally in good agreement except that derived from 
CO, which is outside the 3a range. The high iV(CO)//co ratio gives a upper limit of iV(H 2 ) 
since the CO emission is certainly optically thick. Excluding the value deduced from CO, 
we obtained the average H2 column density of N(fi2) = (7.7 ± 3.8) x 10 21 cm -2 . Hence, the 
error due to the assumptions of structure and optically-thin emission should be less than 
a factor of 2. Combining with the uncertainties caused by various parameters ( excitation 
temperature, distance, source size etc.) and our measurements, we crudely estimate that 
the total error in the calculations of column density and abundance is a factor of ~ 8. The 
column density ratios and the abundance ratios are expected to be more reliable. 
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Results and Discussion 



5.1. Chemistry 

5.1.1. CN and HCN 

Given their relatively large dipole moments, the rotational transitions from the molecules 
CN and HCN are often detected in molecular clouds. HCN can be formed through the disso- 
ciative recombination of HCNH + , while CN is dominantly produced via the photodissociation 
of HCN by UV radiation. Therefore, the CN abundance dramatically increases with in- 
creasing stellar radiation when a star evolves beyond the AGB, through the proto-PN stage 
and to the PN stage. Such a trend has been revealed by a sample study of iBachiller et al. 
(119971 ) . There is observational evidence showing that the abund ance of HCN incre ases in 
the PPN stage, and then dramatically decreases in the PN stage (IHerpin et al.ll2002l ). This 
can be attributed to different roles that dissociative recombination and photodissociation 
play in the two stages. Our observation y ield a HCN/ 13 CO abu ndance ratios of > 0.008, 
consistent the value of > 0.005 reported by IBachiller et al.l (119971 ). However, we do not find 
large CN enhancement in NGC7027. The CN/H 2 abu ndance ratio i s 7.0 x 10~ 8 , about a 
factor of 10 low e r than those in AGB stars reported by I Woods et al.l (120031 ) . The model of 
Hasegawa et al.l (120001 ) seems to overestimate the CN abundance by a factor of 2-3 although 
this may be explained with the uncertainties in the measurements and calculations. The 
unexpectedly low CN abundance was also found by lAli et al.l (120011 ) in their modelling of 
general PNs. 

It has been generally acce p ted t hat the CN/HCN abundance ratio can trace the UV 
radiation field. IBachiller et al.l (119971 ) suggest that after a star evolves into the PN stage, 
the CN/HCN abundance ratio remains constant (~ 9). We obtain iV(CN)/iV(HCN)< 2.4 
in NGC 7027. This probably suggests that the reaction 



CN + H 2 -> HCN + H 



(10) 



is very efficient in this young PN. 



5.1.2. HNC 



HNC can be formed in a similar way as HCN. IHerpin et al.l (120021 ) suggested that 
in PPN stage the HNC abund ance is enhanced with respect to HCN t hrough ion-molecule 
reactions. Both ob servations dSopka et al.lll989l ; IBachiller et al.l Il997l ; iJosselin &: Bachiller 



20031 ) and models (lAli et al.l 120011 ) indicate that the average abundance ratio HNC/HCN 
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in PNs is ~ 0.5. However, we do not detect HNC emission in NGC7027. The HCN/HNC 
ratio is estimated t o be larger than 30. T he ab normally low HNC a b undan ce in NGC 7027 
was also noted by iBachiller et al.l (119971 ) a nd iJosselin fc Bachillerl (120031 ). who obtained 



7V(HCN)/iV(HNC) > 17. As proposed bv IBachiller et al 
destroyed in NGC 7027 through 

HNC + H -> HCN + H, 



(119971 ). HNC may have been 



which is efficient in hi gh temperatures. The h ypothesis is supported by the high HCN abun- 
dance. The model of iHasegawa et al.l (120001 ) did not include the HNC reaction chains. A 
more comprehensive model is required to account for the observed CN/HCN/HNC abun- 
dance ratio. 



5.1.3. HCsN 



HC3N is com m only observe d in AG B stars and proto-PNs (e.g. iCernicharo et al.l 12000 
Pardo et al.ll2004l ). iPardo et al.l (120051 ) found that HC 3 N is quickly reprocessed from HCN 
in PPN stage. Prior to this work, HC3N has never been detected in PNs and its absence is 
often explained as the result o f efficient photodisociation due to the strong UV radiation in 
PNs. iThronson fc Ballv! Jl986h estimated the abundance ratio iV(HC 3 N)/iV(CO) < 1 x 10" 4 
in NGC 7027. The abundance upper limit of HC3N is comparable with that in IRC+10216. 
They thus did not find evidence for underabundance of HC3N in the PN. 

Here we report a possible detection of HC3N in the PN. Our results yield iV(HC3N) /iV(CO) 
1.2 x 10~ 5 , at lest one order of magnitude lower than that in IRC+10216, c onfirming 



the de struction caused by UV photons. The time-dependent chemical model of lAli et al. 



( 120011 ) predicts that the HC3N abundance gradually decreases with the evolution of PNs. 
Their model predicted that NGC 6781, a more evolved PN, has a iV(HC 3 N)/iV(CO) ratio of 
4.2 x 10~ 9 , which is three orders of magnitudes lower than that in NGC 7027. The higher 
abundance of HC3N derived here may be the result of the presence of dense clumps within 
the nebula, which serve to shield the molecules from the UV radiation from the central star. 



5.1.4. C 2 H 



The C2H abundance in NGC 7027 is about a factor of 50 lower than tha t in IRC+10216. 
Based on a comparison between t he observation and the model prediction of IHasegawa et al. 
( 120001 ). IHasegawa &: Kwokl (1200 ll ) concluded that C2H is newly produced via the photodis- 
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sociation of C2H2 instead of being left over from AGB phase. The C2H abundance derived in 
the present study is in good agreement with the value given by their observation and model. 



5.1.5. C 3 H 2 



The C3H 2 abundance in NGC 7027 is at least a factor of 4 lower than that in IRC+10216. 
Albeit with a low abundance, the large dipole mome nt makes C3H9 d etectable. Although 
our detection of C3H2 is uncertain, the observation of ICox et al.l (119871 ) sho wed the presence 



of this species. The abundance of C3H 2 was not predicted by the model of lHasegawa et al. 



( 120001 ) . There is no path for the formation of C3H 2 in PNs, and thus it may be a r emnant 
from AGB envelope. Based on the observation of molecular clouds, iTurnerl ( 119911 ) found 
that CsH 2 avoids the hot environment and favors the cool clouds. Consequently, to explain 
its survival from the dissociation of strong UV radiation field in NGC 7027, we infer that 
C3H2 may arise from some cool high-density clumps. 



5.1.6. HCO+ and N 2 H+ 

HCO + and N 2 H + can be formed through reactions, 

H+ + CO -> HCO+ + H 2 , (12) 

and 

H+ + N 2 -> N 2 H + + H 2 , (13) 

which are expected to have a similar rate coefficient. Both HCO + and N 2 H + can be destroyed 
by dissociative electron recombination. Therefore, if reactions ([121 and (JTBl dominate the 
production of HCO + and N 2 H + , we have the abundance ratio HCO + /N 2 H + m CO/N 2 . How- 
ever, HCO + can be produced through additional processes, such as 

CO+ + H 2 -> HCO+ + H, (14) 

while N 2 H + has no alternative formation channel and can be destroyed at high density by 
proton transfer reaction with CO, 

N 2 H+ + CO -> HCO+ + N 2 . (15) 

As a result, the HCO + /N 2 H + ratio is expected to be larger than the CO/N 2 ratio. 

The observations of HCO + and N 2 H + allow us to test these ion-molecule reactions 
in NGC 7027. Assuming that the N 2 /H 2 abundance ratio is equal to N/H (1.4xl0~ 4 , 
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Zhang et all 120051 ) and N(CO)/N{E 2 ) = 1 x 10" 3 , we have iV(CO)/iV(N 2 ) = 7.1. On 
the other hand, our observations yield iV(HCO + )/iV(N 2 H + ) = 14, higher than the CO/N 2 
abundance ratio by a factor of 2. Hence, additional formation processes for HCO + may be 
significant. 



The N2H + abundance in NGC 7027 is abnormally high, wit h an observed co 



Hasegawa et a" 



umn d ensity 
fl2000h . The 



19971 ). Since 



abundance 



7-8 orders of magnitude higher than predicted by the model of 
species, however, has not been detected in more evolved PNs (IBachiller et al. 
reaction (IT5|) is the unique path for the production of N 2 H + , the high N 2 H H 
therefore implies a high formation rate of Hj~. can be produced through ionization of 
H 2 by cosmic ray or soft X-ray emi ssion from the centra l star. Since strong X-ray emission 
from NGC 7027 has been detected fcastner et alJl200lh . it is possible that X-ray emission 
may play an key role in the formation of N 2 H + in NGC 7027. This is consistent with the fact 
that more evolved PNs generally have weaker X-ray emission and have no N 2 H + emission. If 
high N 2 H + abundance is a property of young PNs, further observations of N 2 H + in a sample 
of young PNs would be useful. 

Although N 2 D + has been extensively observed in dark clouds, we are unable to detect 
it in NGC 7027. From the 3a intensity upper limit of the N 2 D + (l-0) transition, we estimate 
that NGC 7027 has a fractionation ratio iV(N 2 D + )/iV(N 2 H + ) < 0.1. The lower fr actiona- 
tion ratio in NGC 7027 compared to that found in dark clouds (jDaniel et al.l 120071 ) can be 
attributed to different physical environments. 

HCO + in NGC 7027 has a similarly high abundance with that in older PNs. The 
HCO + /HCN abundance ratio in NGC 7027 is . 5, wh ich is in very good agreement with 
the average value for older PNs (IBachiller et al.l 119971 ). The different behaviors of HCO + 
and N 2 H + in young and more evolved PNs suggest that reaction (j!2p is not the dominant 
process for the formation of HCO + . 



5.1.7. HCS+ and CS 



We obtain a iV(HCO + )/iV(HCS + ) column density ratio of 68, which is comparable to 
the ion abundance ratio of N(0 + )/N(S + ) = 111 derived in the ionized regions of NGC 7027 
(IZhang et al.ll2005l ). Hence, the formation rate of HCS + should be similar with that of 
HCO+. 



The estimated CS abundance is about a factor of 10 lower than the model value of 



53 (IHasegawa et al.l 120001 ) . According to iHerpin et al.l (120021 ). after a st ar leaves the AGB 
strong shocks can occur and significantly affect the circumstellar chemistry. IWillacy fc Cherchnefi 
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( 119981 ) suggest that shocks can destroy CS and HCN in circumstellar envelopes, which might 
imply that the low CS abundance observed in NGC 7027 is due to the presence of shocks. 



5.2. Isotopic ratios 

5.2.1. Carbon 

The determination of isotopic abundances plays an important role in our understanding 
of nucleosynthesis in evolved stars. Extensive studies have shown that the 12 C/ 13 C abun- 
dance ratios in red giant star s are considerably lower than the pr ediction made by standard 



stellar evolution models (e.g. Charbonnel fc do Nascimentol fl~998l ) . An extra mixing mecha- 



nism, called cool bottom proce s sing (CBP), was introduced to explain the low 12 C/ 13 C ratio 



pr 
1' 



(ISackmann Sz Boothroydlll999l ; iBoothroyd fe Sackmannlll999l ). For the stars with masses 
in the range 2.5 < M/M & < 6, the hot bottom burni ng (HBB) may oc cur during the AGB 
stage and induce 12 C/ 13 C to further decrease to ~ 3.5 (IFrost et al.lll998T). A few studies have 
been made to study 12 C/ 13 C in molecular envelope of PNs (e.g. iPalla et al.ll2000l ; iBalser et al. 
20021 ) and lend support to the existence of CBP. 



Three 13 C-bearing molecular species, 13 CO, H 13 CN, and H 13 CO + , have been detected in 
our observations, allowing us to estimate the 12 C/ 13 C ratio. The integrated intensity ratio of 
the CO(2-l) and 13 CO(2-l) is 25, yielding an abundance ratio of iV( 12 CO)/iV( 13 CO) = 73. 
Varying values o f the * 2 CO/ 13 CO ratio in NGC 7027 have been re ported in the litera ture. 
Knapp fc Cha^gl fjl985h obtained 40 < iV( 12 CO)/iV( 13 CO) < 130. ISopka et all Jl989h did 
not detect the 13 CO emis si on in NGC 7027 and estima ted 12 CO/ 13 CO to b e large r than 



100. iKahane et all (11992f ). iBachiller et all (I1997f ). and Josselin k Bachiller! (120031 ) got a 



lower limit of 65, 25, and 11, respectively. More recently, Balser et al. ( 20021 ) obtained that 
iV( 12 CO)/iV( 13 CO) ~ 31, about half of our determination. 

The HCN/H 13 CN abundance ratio is 19, similar to the J(HCN(3-2))//(H 13 CN(3-2)) 
integrated intensity ratio of 15. The value may bear a larger error due to the weakness 
of the H 13 CN(3-2) transition. This is for the first time that the isotopomer ratio of this 
species in NGC 7027 is obtained. The value is a factor of ~ 4 lower than the 12 CO/ 13 CO 
ratio. The higher 12 CO/ 13 CO ratio suggests that 13 CO migh t be partly affected by s elective 
photodissociation of the UV radiation field, as proposed by iKnapp fc Changl (119851 ) . Since 
both the CO and HCN lines are likely to be optically thick, the above values may represent 
lower limits. If the 12 C/ 13 C ratio is indeed as low as suggested above, this may imply a 
non-standard mixing process such as CBP. 



Alternatively, the HCO + /H 13 CO + ratio can be used. The HCO + lines are weaker but are 
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definitely optically thin. We used the HCO + (1-0) line at 89.188 GHz and the H 13 CO+ (3-2) 
line at 260.255 GHz since the H 13 CO+(1-0) transition at 88.632 GHz falsi within a spectral 
region which has a poor signal-to-noise ratio and the HCO + (3-2) emission at 267.558 GHz 
falls out of our frequency ranges. Our calculations yield a rather high HCO + /H 13 CO + 
abundance ratio of 117 (see Table [5]). If the HCO + (1-0) transition is optical thick, the 
value is even higher. On the other hand, based on the 3cr upper limit of the integrated 
intensity of the H 13 CO+(1-0) line, we estimate that iV(HCO+)/iV(H 13 CO+) > 60. The 
HCO + /H 13 CO + ra t io der ived in the present study is higher than the value of 40 given by 



Hasegawa &: Kwokl (120011 ). The abnormal carbon isotopic ratio derived from HCO + probably 



suggests that chemical fractionation is significant for this species. 



5.2.2. Oxygen 



During the AGB stage, the nucleosynthesis and dredge-up proc esses cause d estruction of 
18 and enhancement of 17 at the the surface of evolved stars (see Busso 2006, for a recent 



review). iKnapp &: Chang! (119851 ) detected the oxygen isotopic ratio for a sample of carbon- 
rich circumstellar envelopes and found the range of 300-800 and 300-1300 for 16 0/ 17 and 
16 0/ 18 0, respectively. Based on the fact that the 17 0/ 18 ratios in carbon-rich envelopes 
are higher than that in the interstellar medium by a factor of 4-5, they ruled out AGB stars 
as the dominant source of oxygen in the interstellar medium. 

From the detection of CO, C 17 0, and C 18 emission in NGC7027, we obtain the oxy- 
gen isotopic ratios iV( 16 0)/iV( 17 0) = 698, iV( 16 0)/iV( 18 0) = 1763, and iV( 1 7 OW( 18 0) = 3 
Our r esults are consistent with the ratios (or the lower limits) derived by IKnapp fc Chang 
Jl985h in NGC 7027. Compared with the oxygen isotopic r atios in the Sun ( N( 16 0) /N( 17 0) = 
0.2, iV( 16 0)/iV( 18 0) = 499, and iV( 17 0)/iV( 18 0) = 2682. lLoddersll2003h . NGC 7027 clearly 
shows enhancement of 17 and destruction of 18 O. which well agrees with the predictions of 
stellar models. 



5.3. Hydrogen and helium recombination lines 

A total of over 30 H and He recombination lines are detected in our observations. These 
include nine Una lines, nine Hn/3 lines, and probably six B.n'-f lines and eight Hena lines. 
The linewidths and integrated intensities of these lines are listed in Table HJ 



Due to optical depth effects, the line widths and radial velocit ies of recombination 
lines at different frequencies may be different (jErshov fc Berulis Ill989l ). The high frequency 
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emission has negligible optical depth, and is likely to be emitted from the whole ionized 
region. The low frequency emission, however, is more likely to be optically thick and if so will 
be emitted from the outer regions. Our data do not show any clear correlation between the 
line frequencies and the line widths and velocities, suggesting that these recombination lines 
are optically thin. From profile fittings, we determine the averaged line widths (FWHM) for 
Una, Hn/3 and Hena lines to be 44±8kms _1 , 54±9kms _1 , and 38±9kms _1 , respectively. 
No systematic difference is found for the widths of hydrogen lines and the helium lines. The 
Una lines detect ed here are a factor o f abo ut 1.6 narrower than the low frequency HllOa 
line detected by IChaisson fc Malkan Jl976h . which is expected to be optical thick. The 
observed widths o f hydr ogen radio recombination lines can be interpreted by the model of 
Ershov fc Berulis I Jl989h (cf. their Fig.2c). 



The relative intensities of observed Una lines as a function of their quantum numbers 
are plotted in Fig. [81 In this plot, we have corrected for the effect of beam dilution and 
translated the inten s ity un it into Jy. The low frequency lines H76a and HI 10a reported by 
Chaisson fc Malkan! (119761 ) are also plotted in this figure. We have normalized the intensities 
of the H76a and HllOa lines to J(H39a;) = 1 by assuming that they have the same beam- 
filling factor, which may intr oduce an error less than a factor of 5. For comparison, we plot 
the theoretical predictions of IStorey fc Hummer Jl995h under the assump tion of Case B with 
severa l combinations of electron temperatures and densities. Note that IStorey &: Hummer 
( 119951 ) only gives the line emissivities for n < 50 and the predictions for n > 50 in Fig. [8] 
are linear extrapolations of their results. If a correction for continuum free-free opacity 
is applied, the predicted lines would have a larger slope than those shown in Fig. [SJ An 
inspection of Fig. [8] shows that the intensities of the high-frequency Una lines (> 50 GHz) 
are consistent with the theory of recombination and are not affected by free-free opacity. 
The fact that the high frequency Una lines are weaker than the theoretical predictions may 
suggest that they are optically thick. 

Una and Hn0 radio recombin ation lines were also detected in molecular clouds (e.g. 
Turner! 1 1 99 ll ; iNummelin et al.ll2000i ). allowing us to make a comparison study. The averaged 
main-beam temperatures for Una lines and Hn/3 lines determined from our observations are 
T(a) = 0.046 K and T(j3) = 0.014 K, respectively. The corresponding T R (a)/T R (j3) ratio 
in NGC7027 is 3.3, which is in g ood agreement with that fou nd in two molecular clouds 
OMC-1 and SgrB2 ([Turner! 1 19911 ). According to iTurneri (119911 ) . stimulated emission must 
be significant in order to explain such a Tr(«) /Tr(/3) ratio. The observed T R (a) value in 
NGC 7027 is lower than that in OMC-1 by a factor of ten, whic h can be attrib uted to the 
smaller size (~ 0.1 pc) of this PN. Following the calculations of ITurneri (Il99ll ). NGC 7027 
has an emission measure EM ~ 6 x 10 6 cm -6 pc. 
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Thum et al. I (119981 ) studied the radio recombination lines of hydrogen in the emission- 
line star MWC 349 and found that the amplification factor of the recombination line maser 
has a peak value near n = 19. Fig. [S] compares hydrogen recombination lines in MWC 349 
and those in NGC7027. We can clearly see that for MWC 349, the intensities of the Una 
lines near n = 30 are consistently higher than the theoretical pred ictions, which has been 
attributed to the broadband maser phenomenon (IThum et al. Ill998l ) . The trend is not found 
in NGC7027. From the T{a)/T{(3) ratio, we find stimulated emission to be significant for 
the hydrogen radio recombination in NGC 7027 and infer that all the Una lines observed in 
the present work have a similar amplification factor. NGC 7027 has a density of 10 4 cm~ 3 , 
ab out a factor of 10 4 lower than that of MWC 349. Consequently, according to the model 
of IThum et al. I (119981 ). the quantum number of the peak amplification (n max ) in NGC 7027 
is larger than that in MWC 349. Since n max = 19 ± 2 in MWC 349, a similar amplification 
factor for the Una would imply that n max is in the range between 34-45 in NGC 7027. 

The hydrogen and heli um radio recombina tion lines can also be used to calculate the 
He + /H + abundance ratio. iPalmer et al.l (119691 ) deduced the He/H abundance ratio from 
radio recombination for a sample of H II regions and f ound that the resul ts are in good 
agreement with those derived from optical data. Following IPalmer et al.l (119691 ) and assuming 
that the ratio of the integrated intensity of the Hena line to that of the Una represents the 
He + /H + abundance ratio, we obtain iV(He + )/iV(H + ) = 0.11 ±0.02. An additional 20% error 
is given by uncertainties in the measurements. The derived He + /H + abund ance ratio is 



abou t a factor of two higher than that determined from the optical spectra (jZhang et al. 



20051 ). A possible explanation is that the He + zone is more extended than the H + zone so 
that the effect of beam dilution is more severe for the H I radio recombination lines. To 
account for the discrepancy between the abundance derived from the radio data and that 
from the optical data, the size of the He + zone may be larger than the H + zone by a factor of 
about 1.3. This is possible if the ionization spectrum is hard enough and more UV photons 
have an energy close to the ionization potential of helium. In that case, the He 2+ zone is 
expected to be large in the high excitation PNs. However, the He II radio recombination lines 
are below our detection limit, and thus we have not attempted to calculate the elemental 
abundance ratio of helium to hydrogen. 



5.4. Line profiles 



Fig. [5] and Fig. [6] present the line profiles detected with the ARO 12m telescope and 
the SMT 10m telescope, respectively. The HCO + line exhibits a self-absorption featur e 
at Vlsr ~ 22km/s, in accord with the observations of HCO + line by iDeguchi et al.l (Il990l ). 
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We also find a weak absorption for the HCN lines at Vlsr ~ 30km/s. 

Table. [T] gives the FWHM of detected molecula r lines. The emission lin es from CO 
and its isotopomers have a width between 23-30 km/s. lHasegawa & Kwokl (120011 ) found that 
the HCN, CN, and C2H have larger line widths than that of the CO emission, even after 
taking into account the effects of hyperfine components. Their findings are confirmed by 
our observations. One possible explanation for the larger linewidths is that these molecular 
lines contai n components arising from fast-moving clump s within the neutral envelope of 
NGC7027 (IRedman et al. 1120031 ; iHuggins fc Mauronll2002l ). The CO emission originates in 
a much more extended region compared to other molecular species, and thus is less affected 
by the contribution from the fast components. The present observations show that the 
lines from HCO + and N 2 H + are consistently narrower than that from CO. We therefore 
infer that the fast clumps might be neutral and have no contribution to the emission from 
the molecular ions. The molecular ions mostly originate from the PDR, which has a lower 
expansion velocity. Further observations of lines from molecular ions are required to verify 
this point. Fig. [5] and Fig. [6] show that the recombination lines are generally broadener than 
CO lines, suggesting that the acceleration by the UV radiation field is significant in the 
ionized regions of NGC 7027. 



5.5. Comparison with the spectra of IRC+10216 and CRL 618 



Comparison with the spectra of AGB stars, PPNs, and PNs can provide essential infor- 
mation of chemical evolution in la te-type stars. Rec e nt mo lecul ar line survey s of IR C+10216 
and CRL 618 were presented by ICernicharo et al.l (120001 ) and iPardo et al.l (120071 ). respec- 
tively. The frequency ranges are 129-172.5 GHz for IRC+10216, and 80.25-115.75 GHz, 
131.25-179.25 GHz, and 204.25-275.25 GHz for CRL 618. Since these spectra were obtained 
using different instruments, it is hard to make a quantitative comparison with our observa- 
tion of NGC 7027. Through a qualitative comparison, we summarize the main differences 
between the spectra of the three objects as follows: 



1. The molecular lines detected in NGC 7027 are much fewer and generally fainter than 
those in IRC+10216 and CRL 618, even after the correction of the relative distances of 
the three objects. This can be attributed to a) the low density in the neutral envelope 
of NGC 7027, whic h makes it harder to colli sionally excite the higher rotational levels 
of most molecules (IThronson &: Ballyl Il986l ). and b) destruction of molecules caused 
by the strong UV radiation and shock waves in NGC 7027. 



2. Refractory Metal-bearing and silicon-bearing species which are plentifully present in 
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IRC+10216 have only week emission in CRL618 and are not detected in NGC7027. 
These molecules may be depleted onto dust grains with stellar evolution. 

3. The CS emission is strong in IRC+10216, medium in CRL618, and not seen in 
NGC 7027. This might suggest that shock waves play an important role in the chem- 
istry of NGC 7027. 

4. Ionized species (HCO + , HCS + , and N2H + ) and recombination lines detected in CRL 618 
and NGC 7027 are not seen in the spectra of IRC+10216. This is due to photoionization 
by the hotter central stars of CRL 618 and NGC 7027. 

5. The N 2 H + /HCO + abundance ratio in NGC 7027 is about a factor of 30 higher than 
that in CRL 618, supporting the hypothesis that H3 is dominantly produced by soft 
X-ray from the central star. 

The comparison between the spectra of CRL 618 and NGC 7027 is particularly mean- 
ingful because the former is a much evolved PPN and will evolve to become a young PN 
like NGC 7027 after only a few hundred years. From the above comparison, we can conclude 
that the chemical compositions can dramatically change after the PPN stage. 

6. Summary 

In this spectral-line survey of the PN NGC 7027, we have detected a total of 67 spectral 
lines consisting of 21 molecular lines, 32 hydrogen and helium recombination lines, and 14 
unidentified line. The line intensities were used to calculate excitation temperatures, column 
densities, and fractional abundances of the observed molecular species. The main findings 
are summarized as follows: 

1. The low CN/HCN ratio seems to go against the theoretical expected photodissociation 
of HCN into CN in a strong UV environment. This could imply that the photodisso- 
ciation process might be compensated by the unexpectedly efficient reaction between 
CN and H2. Alternatively, a large amount of HNC may have been transferred into 
HCN in the hot envelope of NGC 7027. 

2. From the HCO + /N 2 H + abundance ratio, we suggest that the reaction of CO + and H 2 
may be important for the production of HCO + . The soft X-ray emission from the 
central star play to important role for the production of N 2 H + in NGC 7027. 
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3. We find evidence for underabundance of CS, suggesting that the species may have been 
destroyed by shocks in the young PN. 

4. Although the estimated CO/ 13 CO, HCN/H 13 CN, and HCO+/ 13 CO+ ratios are differ- 
ent, the resulting low 12 C/ 13 C ratio supports the presence of extra mixing process in 
AGB stage. The oxygen isotopomer ratios is consistent with the predictions of stellar 
models. 

5. The intensities of hydrogen recombination lines suggest that stimulated excitation may 
be significant. Our estimates of He + /H + abundance ratio from the ratio recombination 
lines is higher than that from the optical lines. The discrepancy can be explained if 
the He + zone is more extended that the H + zone in the young PN. 

6. The different widths of lines from CO, other neutral molecules, molecular ions, and 
hydrogen and helium can be used as dynamical probes of the circumstellar envelope of 
NGC 7027. Our results suggest the presence of fast neutral clumps within the molecular 
envelope of NGC 7027. 

These observations show that the circumstellar envelope of NGC 7027 can serve as a very 
useful chemical laboratory. Comparison with existing theoretical chemical models suggests 
that a more comprehensive model is required to explain the observed molecular abundances, 
in particular for the production of N2H + , the destruction of CS, and the survival of HC3N 
and C3H2 in the PN. A more comprehensive model will need to include effects of soft X-ray 
emission from the central star, shocks, and dense clumps. 

Changing physical conditions are the main factor for the different molecular abundances 
in the circumstellar envelopes of stars at different evolutionary stages. The observational 
results reported in the present paper provide a foundation for further studies of circumstellar 
chemistry. In this spectral survey program, we also obtained radio spectra of two AGB stars, 
IRC+10216 and CIT6, and one proto-PN, CRL2688, in addition to the PN NGC 7027. A 
systematic comparison of the chemical compositions in different evolutionary stages will be 
presented in a forthcoming paper. 

We thank Aldo Apponi for his help in the processing of the ARO data. We also thank 
Jun-ichi Nakashima and Jin-Hua He for useful discussions. This work is supported in part by 
a grant awarded to SK from the Research Grants Council of Hong Kong. DVT acknowledges 
support of this work by Academia Sinica, Taiwan, and the Natural Science Council of Taiwan. 
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Table 1. Detected molecular transitions. 



Species 


Transition 


Frequency 


rms 


R 


/ T n dv a 


FWHM a 


Remarks 6 




(upper-lower) 


(GHz) 


(mK) 


(K) 


(K km/s) 


(km/s) 




CO 


2-1 


230.538 


3.0 


5.824 


141.97 


24.2 


H01,T86,B97 


13 CO 


1-0 


110.201 


3.3 


0.069 


1.95 


27.5 


S89,B97 




2-1 


220.398 


3.7 


0.204 


5.72 


29.1 


H01,J03,B97 


C ls O 


2-1 


219.560 


3.6 


0.020 


0.41 


23.3 




c 17 o 


2-1 


224.714 


4.3 


0.032 


0.86 


25.4 




CN 


2-1,3/2-3/2 


226.314 


3.9 


0.028 


1.79 




*,C93,J03,B97 




2-1,3/2-1/2 


226.659 


3.9 


0.091 


4.95 


44.5 


*,C93,J03,B97 




2-1,5/2-3/2 


226.874 


3.9 


0.197 


7.76 


36.5 


*,C93,J03,B97 




2-1,5/2-1/2 


227.192 


3.3 


0.008: 


0.26: c 






C 2 H 


3-2,7/2-5/2 


262.004 


5.5 


0.041 


2.01 


94.1 


*,H01 




3-2,5/2-3/2 


262.067 


5.5 


0.033 


1.74 


96.7 


*,H01 


C3H2 


432-321 


227.169 


3.3 


0.008: 


0.26: c 








625-5l4 


251.527 


4.0 


0.010: 


0.19: 








533-422 


254.988 


4.1 


0.012: 


0.20: 






HCN 


1-0 


88.632 


3.8 


0.059 


1.94 


36.6 


D86,D90,S89,C93,J03,B97 




3-2 


265.886 


5.5 


0.207 


8.22 


56.6 


H01,C93 


H 13 CN 


1-0 


86.339 


3.5 


0.009: 


0.28: 








3-2 


259.012 


4.3 


0.013 


0.56 


64.9: 




HCO+ 


1-0 


89.188 


3.7 


0.247 


5.39 


19.2 


S89,D90,C93,J03,B97 


H 13 CO+ 


3-2 


260.255 


4.1 


0.017 


0.27 


16.1 


H01 


HC 17 0+ 


1-0 


87.057 


3.4 


0.008: 


0.21: 






HC 18 0+ 


1-0 


85.162 


7.4 


0.015: 


0.30: 






HCS+ 


6-5 


256.028 


4.0 


0.014 


0.25 


17.0: 


H01 


HC3N 


8-7 


72.783 


3.6 


0.014 


0.34 








9-8 


81.881 


3.6 


0.015 


0.31 








10-9 


90.979 


3.1 


0.008: 


0.16: 








11-10 


100.076 


2.9 


0.011 


0.22 








12-11 


109.174 


2.8 


0.004: 


0.12: 








26-25 


236.513 


3.4 


0.004: 


0.11: 








27-26 


245.606 


2.7 


0.014 


0.23 








28-27 


254.700 


4.1 


0.011: 


0.13: 






N 2 H+ 


1-0 


93.173 


2.6 


0.017 


0.35 


19.7 


C93,J03 



a the symbol " :" indicates uncertain detections. 

b Also detected by: B97-Bachiller et al. (1997); C93-Cox et al. (1993); D86-Deguchi et al. (1986); D90-Deguchi et 
al. (1990); J03-Jossclin & Bachiller (2003); HOl-Hascgawa & Kowk (2001); S89-Sopka ct al. (1989); T86-Thronson & 
Bally (1986). *-unsolved fine-structure lines. 



c The CN at 227.192 GHz and the C 3 H 2 at 227.169 GHz lines are blended with each other. 
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Table 2. Unidentified features. 



Frequency rms 


T R 


jT R dv 


(GHz) (mK) 


(K) 


(K km/s) 



74.739 


3.8 





012 


0.63 


89.656" 


3.6 





012 


0.26 


90.739" 


3.1 





013 


0.23 


95.246 


3.2 





Oil 


0.47 


96.236" 


5.3 





019 


0.35 


98.079 


2.9 





010 


0.20 


102.789 


3.3 





013 


0.29 


107.802 


2.8 





012 


0.37 


229.256 fe 


3.7 





015 


0.22 


257.637 6 


4.8 





015 


0.25 


261.569 6 


5.5 





020 


0.28 


261.820 


5.5 





019 


0.36 


263.369 6 


5.0 





017 


0.26 


264.744 


5.2 





019 


0.36 



a Also detec ted in Sagittarius B2 (M) 
humeri Il989l ): 



b Also detectec 


in S< 


(Nummelin et al. 


1998) 
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Table 3. Non-detected molecular species in NGC 7027. 



Species 


Transition 


Frequency 


rms 


N( 


:m 2 ) 


AT(obs.)/ 




(upper-lower) 


(GHz) 


(mK) 


obs." 


modcl b 


N (model) a 


CS 


v=Q J=5-4 


244.936 


2.7 


1.9cl2 


2.5cl3 


0.08 


C 3 H 


2 n 3/2 J=23/2-21/2 b 


263.332 


4.9 


l.lel3 






C 4 H 


N=24-23 a 


228.349 


2.9 


1.7cl3 


1.4ell 


121 


SiC 


3 n 2 J=6-5 


236.288 


3.5 


7.5cl2 


2.6cl2 


2.9 


SiO 


u=0 J=6-5 


260.518 


5.2 


1.0cl3 


1.5el3 


0.7 


SiS 


u=0 J=14-13 


254.103 


3.7 


8.9cl3 






HNC 


J=l-0 


90.664 


3.4 


1.4el2 






HC 5 N 


J=27-26 


71.890 


3.6 


1.5el2 






CH 3 CN 


J K =5(1)-4(1),5(0)-4(0) 


91.986 


2.4 


2.3el3 






H 2 CO 


3(1,2)-2(1,1) 


225.678 


3.9 


4.5el2 


9.2e9 


489 


N 2 D+ 


J=l-0 


77.109 


4.1 


2.8cl2 







a Upper limits. 

b Calculated by J n(i)4:-Kr 2 dr/AS from Hasegawa ct al. (2000), where AS is the projected 
source area. 
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Table 4. Recombination lines in NGC 7027. 



Line 


Frequency 


Tr 


/ T R dv 


FWHM 




(GHz) 


(K) 


(K km/s) 


(km/s) 


H44a 


74.645 


0.045 


2.01 


46.6 


H43a 


79.913 


0.037 


1.77 


52.5 


H42a 


85.688 


0.037 


1.68 


35.3 


H41a 


92.034 


0.054 


2.39 


41.2 


H40« 


99.022 


0.051 


2.19 


43.5 


H39a 


106.737 


0.049 


2.47 


49.4 




1 fin 919 




9 77 


AO fi 


H30a 


231.901 


0.037 


1.88 


46.2 


H29« 


256.302 


0.053 


2.46 


44.3 


H55/3 


74.940 


0.013 


0.46 


52.9 


H54/3 


79.104 


0.025 


0.68 


35.7 


H53/3 


83.582 


0.011 


0.41 


48.5 


H52/3 


88.406 


0.016 


0.60 


60.1 


H51/3 


93.607 


0.010 


0.40 


64.5 


H50/3 


99.225 


0.013 


0.53 


49.8 


H49/3 


105.302 


0.012 


0.57 


65.5 


H37/3 


240.021 


0.013 


0.54 


57.5 


H36/3 


260.033 


0.011 


0.42 


92.1: 


H6I7 


80.900 


0.012 


0.48 


50.6 


H6O7 


84.914 


0.024 


0.29: 


17.8: 


H587 


93.776 


0.009 


0.28 


64.1 


H57 7 


98.671 


0.006 


0.16: 


31.0: 


H567 


103.915 


0.009 


0.22: 


35.5: 


H557 


109.536 


0.008 


0.25 


35.1 


He44a 


74.674 


0.009 


0.21 


40.2 


He43a 


79.945 


0.012 


0.17: 


33.3: 


He42a 


85.723 


0.023 


0.37: 


36.3: 


He41a 


92.072 


0.012 


0.28 


27.0 


He40o 


99.063 


0.008 


0.22 


55.5: 


He39o 


106.781 


0.005 


0.18: 


31.6: 


He30o 


231.995 


0.009 


0.23 


42.0 


He29o 


256.406 


0.014 


0.35 


35.2 
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Table 5. Excitation temperatures (T cx ), column densities (N), and abundances relative to 

H 2 (fx). 





T (TO 


iV(cm- 2 ) 


TV (obs 1 / 


JX 


N/fv(cm" 2N ) 






obs. 


model a 


N(xr\nc\p\\ 

1 V I 111UUC1 1 






HC 3 N 


34.6 








8.5e-9 


5.2e21 


CN 




1.02el4 


2.4el4 


0.43 


7.0e-8 


1.5e21 


C 2 H 




3.49el4 


1.5el4 


2.33 


5.4e-8 


6.5e21 


C3H2 




l.llel4: 






8.3e-9: 


1.3e22: 


HCS+ 




4.96el2 


9.0ell 


5.51 


1.0e-9 


5.0e21 


N 2 H+ 




2.43el3 


3.2e5 


7.59e7 


3.8e-9 


6.4e21 


HCO+ 




3.38el4 


4.6el3 


7.35 


4.8e-8 


7.0e21 


H 13 CO+ 




7.23ell 






4.1e-10 


1.8e22 


HC 17 0+ 




2.85el3: 






2.3e-9: 


1.2e22: 


HC 18 0+ 




3.89el3: 






3.2e-9: 


1.2e22: 


HCN 


5.4 


7.05el4 


1.8el4 


3.92 


4.5e-8 


1.6e22 


H 13 CN 




3.65el3 






3.9e-9 


9.4e21 


CO 


40.1 


3.79el8 


3.2e20 


0.01 


l.le-4 


3.4e22 


13 CO 




5.19el6 






7.4e-6 


7.0e21 


c 18 o 




2.15el5 






3.8e-7 


5.7e21 


c 17 o 




5.43el5 






7.4e-7 


7.3e21 



a Calculated by J n(i)A7rr 2 dr / AS from Hasegawa et al. (2000), where AS 1 is 
the projected source area. 

b For the species with optically thick emission, this is a lower limit (see text 
for details). 
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Fig. 1— From top to bottom: the 71-111 GHz 157-161 GHz and the 157-161 GHz MAC 
spectra of NGC 7027 obtained with the ARO 12 m telescope. The spectra have been 
smoothed to a resolution of 1 MHz. The spectra at full resolution is available online in 
Fig. El 




Fig. 2.— The 218-267 GHz FFB spectrum of NGC7027 obtained with the SMT 10 m tele- 
scope. The spectra have been smoothed to a resolution of 3 MHz. The spectra at full 
resolution is available online in Fig. HI 




Fig. 3.— The 71-111 GHz and 157-161 GHz spectra of NGC7027 obtained with the ARO 
12 m telescope. The black and red lines are the MAC and FB data, respectively. The spectral 
resolution is 1MHz. The ':' represents uncertain detection. 
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Fig. 4.— The 218-267 GHz spectrum of NGC7027 obtained with the SMT 10 m telescope. 
The black and red lines are the FFB and AOS data, respectively. The spectral resolution is 
3 MHz. The ':' represents uncertain detection. 
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Fig. 5. — Observed line profiles of NGC7027 by the ARO f2m telescope. Positions and 
relative intensities of hyperfine components are marked by vertical solid lines. The spectral 
resolution is 1 MHz. The velocity range of 13 CO emission (above la) is indicated by vertical 
dashed lines. 
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Fig. 5. — Continued. 
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Fig. 6. — Observed line profiles of NGC7027 by the SMT 10 m telescope. The description 
is same with those in Fig. [5] The spectral resolution is 3 MHz. 
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Fig. 6. — Continued. 
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Fig. 7. — Rotational diagram for HC 3 N in NGC 7027. The solid line is obtained through a 
linear least-squares fit of all the data. 
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Fig. 8. — Relative fluxes of H I recombination a-transitions in MWC 349 (open circles; 
Thum et al. 1998) and NGC 7027. The open triangles and the filled triangles are from 
our observations of NGC 7027 with the ARO 12m and the SMT 10m telescopes, respectively. 
The error bar on the lower left indicates typical uncertainties of our obervations. The 
filled circles are from the observations of NGC 7027 by Chaison & Malkan (1976). The 
theoretical predictions of Storey and Hummer (1995) are given by assuming T e = 12500 K and 
N e = 10 4 cm" 3 (solid line), T c = 12500 K and N c = 10 3 cm- 3 (dotted line), and T e = 10000 K 
and N e = 10 4 cm -3 (dashed line). No correction for free-free continuum opacity is included. 



